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Abstract: We consider the application of quantum corrections computed using renormal- 
ization group arguments in the astrophysical domain and show that, for the most natural 
interpretation of the renormalization group scale parameter, a gravitational coupling pa- 
rameter G varying 10 -7 of its value across a galaxy (which is roughly a variation of 10~ 12 
per light-year) is sufficient to generate galaxy rotation curves in agreement with the obser- 
vations. The quality of the resulting fit is similar to the Isothermal profile quality once both 
the shape of the rotation curve and the mass-to-light ratios are considered for evaluation. 
In order to perform the analysis, we use recent high quality data from nine regular disk 
galaxies. For the sake of comparison, the same set of data is modeled also for the Modified 
Newtonian Dynamics (MOND) and for the recently proposed Scalar Tensor Vector Gravity 
(STVG). At face value, the model based on quantum corrections clearly leads to better fits 
than these two alternative theories. 

KEYWORDS: dark matter theory, rotation curves of galaxies, quantum field theory on 
curved space. 



Contents 



1. Introduction 



1 



2. Effective action of gravity and its low-energy behaviour 

2.1 Covariance arguments and the scale setting 

2.2 The dynamics for a slowly varying gravitational coupling parameter 



4 



4 



7 



3. Galaxy rotation curves 

3.1 General properties 

3.2 Mass decomposition: our results on real galaxies in comparison with other 
models 



11 



11 



13 



4. Conclusions and discussions 



19 



1. Introduction 

The dynamics of many galaxies looks like as if the main part of their masses is distributed 
in a different way than their light is, with heavier densities than those expected from 
the emitted radiation by the stars and gas. Curiously, although galaxies are far from 
being "hairless" , like fundamental particles or black holes, these masses discrepancies follow 
considerably strong patterns, including a series of correlations with the luminous matter, 



Usually, in order to solve this missing matter problem, a qualitatively new kind of 
matter, dark matter (DM), is evoked. Among all kinds of dark matter, current cosmologi- 
cal observations (like the large scale structure and the microwave background anisotropies) 
favor the collisionless and "cold" type of dark matter (CDM). However the CDM paradigm 
faces some difficulties at galactic scales, like the angular momentum issues on galaxy for- 
mation, the discrepancies on the number of satellite galaxies, and the cuspy dark matter 
density profile (see [3] for a recent review). The CDM paradigm has not yet yielded a 
satisfactory dark matter profile for spiral galaxies from the simulations of cosmological 
evolution. The problem with the cuspy dark matter distribution in galaxies was already 
pointed in [4, 5], and in particular both the Navarro- Frenk- White [6] and the Moore [7] dark 
matter profiles have a cusp on the dark matter density at the galactic center whose dynam- 
ical consequences do not match the observations; the discrepancies with the observations 
becoming particularly clear for the Low Surface Brightness galaxies. The problematics 
associated with this cuspy dark matter profile, and possible solutions, were systematically 
analyzed in many references, including [8, 9, 10, 11, 12, 13], see [14] for a recent review. 



e.g. [1, 2]. 
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It is well known that the parameters of the dark matter profiles that best fit the 
observed rotation curves (like the Isothermal [15] or the Burkert 1 [16] ones, which have 
a core and were proposed on phenomenological grounds) display many correlations with 
the baryonic matter behavior, e.g. [17, 18, 13, 15, 16], see also [19]. These correlations 
between a galaxy rotation curve and its baryonic matter suggest another interpretation 
to the missing mass problem in galaxies, namely a change of the gravitational law itself 
instead of introducing DM. This has lead to diverse proposals, including modifications 
on the law of inertia [20] and its extensions to General Relativity (e.g., [21]), proposals 
with extra-dimensions [22], actions of General Relativity with an extra affine connection 
or a new graviton [23], conformal gravity [24], non-symmetric metrics [25], among others 
possibilities. MOND [20, 26] is by far the most studied example of the last approach. 
Although it seems to be consistent with the general features of the galactic dynamics, its 
concordance with observations (if the original recipe is applied) started to be questioned 
in the light of the higher precision recent observations and the improvements on the stellar 
mass-to-light ratios constraints, e.g. [10, 27, 28, 29, 30]. 

One can suppose that the deviation from the gravitational (either Einstein or Newton) 
law is due to the quantum effects, such as semiclassical corrections, effects of quantum 
gravity, consequences of string theory physics, extra dimensions, branes or some other 
(maybe yet unknown) model of "quantum gravity" . 

The corrections to the Newton law is a relatively common feature of the different models 
of "quantum gravity", including higher derivative quantum gravity [31], different versions 
of the effective low-energy quantum theory of gravitational field [32, 33, 34], so-called non- 
perturbative quantum gravity based on the hypothesis of the existence of the non-Gaussian 
UV fixed point [35] and, finally, in the semiclassical approach to quantum gravity [36] (see 
also references therein) . The application of these corrections has been elaborated recently 
in the cosmological [37, 38, 39] and astrophysical [39] areas. One can see the recent paper 
[40] for the detailed discussion of the quantum field theory backgrounds of these quantum 
effects. Let us emphasize, from the very beginning, that in the most cases the present day 
state of art in all mentioned approaches to quantum gravity does not enable one to really 
calculate the relevant quantum contributions to the Newton law in a unique and consistent 
way. At the same time, we have many reasons to believe that these quantum corrections 
can be non-zero and, in principle, may have a measurable effect. Typically, the theoretical 
estimate for the quantum contribution to the gravitational law involves more or less strong 
arbitrariness, related to the dependence of the quantum corrections on some dimensional 
parameter (scale parameter in the renormalization group-based approaches, for instance) 
from one side, and to the physical interpretation of this parameter from the other side. 

In the present paper we will consider a relatively general model of quantum corrections 
which is mainly controlled by covariance. Furthermore we will use a new identification of 
the scale parameter [i in the astrophysical setting (rotation curves problem), which we 
believe to be the more natural compared to the one (/x ~ 1/r) considered before in 

Considering the CDM simulations, the Burkert profile has merits over the Isothermal one, in particular 
for its large R behavior and its finite total mass; see also [70]. 
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[41, 39, 35] (see also [42]). It turns out that this new identification of scale leads to the 
surprisingly good result for the rotation curves of the galaxies, which is very close in quality 
to the output of the mainstream approach based on the DM paradigm. 

The previous papers on the subject were based on the renormalization group equa- 
tions coming from the higher-derivative quantum gravity [41] (unfortunately, the proper 
renormalization group equations [31] are ambiguous in this case, making their applica- 
tion to astrophysics doubtful [37]), on the assumption of the Appelquist and Carazzone 
decoupling of the quantum corrections in the low-energy domain [43, 37, 38] and on the 
assumption of the asymptotic safety and existence of the non-Gaussian fixed point in the 
four-dimensional quantum gravity [35]. It is remarkable that all those, in fact rather dif- 
ferent approaches, converge in predicting qualitatively similar logarithmic running of the 
effective Newton constant with the scale \i. It was shown in [39] (see also previous qual- 
itatively similar consideration in [41]) that this logarithmic behaviour, together with the 
mentioned above identification of pT x to the distance from the center of the galaxy, r, can 
partially explain the main features of the rotation curves without invoking the DM. The 
astrophysical applications of [41], [39] and [35] were performed for the point-like model of 
the galaxy, which is not supposed to be realistic. It was pointed out in [39] that it would 
be very interesting to make a more detailed investigation, taking the case of a plane or thin 
disk distribution of mass in the galaxy. We present the corresponding analysis and show 
that it meets certain obstacles at both theoretical and phenomenological levels. 

In order to make a conclusive consideration of the astrophysical application of the 
renormalization group method, one has to ask whether the /i. ~ \jr identification is the 
unique possible one, or there are some alternative options. One can remember that the 
identification of scale in the cosmological setting is a nontrivial issue, which was treated 
differently by different authors in different papers [44, 37, 43, 38, 45, 35]. Furthermore, 
was an interesting attempt to construct the regular scale- fitting procedure [46]. The output 
of this procedure was the most natural identification of \i with the Hubble parameter H, 
which is the energy characteristic of the external metric leg of the Feynman diagram, 
corresponding to the quantum correction to the expansion of the universe. Indeed, similar 
energy characteristic can be constructed in the case of the galaxies rotation curves, and it 
is indeed different from 1/r. Moreover, this "right" setting of the energy scale leads to the 
immediate and dramatic improvement of the phenomenological analysis results. After all, 
we gain a chance to explain the rotation curves on the basis of quantum corrections to the 
classical action of gravity, without using DM. 

The paper is organized as follows. In section 2 we discuss very general features of quan- 
tum contributions and establish their likely form, which is based on the general covariance 
and on the proper existence of these contributions, following [39]. As we have already 
mentioned above, this existence can not be either proved or disproved at the present-day 
state of art in this area [40] and it is legitimate to apply a phenomenological approach. Let 
us note that the covariance enables one to put serious restrictions on the possible forms 
of quantum corrections in the cosmological setting, the astrophysical ones can be deduced 
starting from the cosmological setting result and the assumption of a unique effective ac- 
tion for gravity. The main new element in section 2 is the new definition of the energy 
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scale, which is appropriate for the application to the rotation curves problem. In section 3 
we present the numerical results of our model applied to the data of nine regular galaxies 
divided into two samples of data. For comparison purposes, our results are shown together 
with the numerical results of three other proposals, obtained from exactly the same data 
and procedures. These three proposals are: dark matter (modeled with the Isothermal pro- 
file [15]), MOND [20] and STVG [56, 25]. The first is the leading candidate for the missing 
matter problem in galaxies and all the universe altogether, MOND is the paradigmatic 
model for galaxy rotation curves without dark matter and STVG is a recently proposed 
model for astrophysics and cosmology without dark matter that also uses the variation of 
the gravitational coupling parameter G to fit galaxy rotation curves. Finally, in section 4 
we draw our Conclusions and discuss the perspectives of the theory of quantum corrections 
(especially the ones based on the renormalization group) as a potential competitor of the 
standard ACDM model, at both astrophysical and cosmological scales. 

2. Effective action of gravity and its low-energy behaviour 

In this section we will discuss the possible form of Effective Action of gravity (EA), including 
the contributions coming from quantum matter or from the quantization of the proper 
gravitational field. The introduction to the gravitational theory at quantum level can be 
found in [47, 36] (see also recent reviews in [48] and [40]). 

2.1 Covariance arguments and the scale setting 

The EA is a functional which is the analog of the classical gravitational action at quantum 
level. The EA can be seen as a classical action plus quantum contributions. The finite 
part of EA is a non-local functional, which usually can not be calculated explicitly. Let us 
use some general features of EA to establish a possible form of low-energy quantum effects 
at the cosmological and astrophysical scales [37, 38, 39, 48]. Our assumptions include the 
covariance of EA, and that the low-energy gravity should not have other light degrees of 
freedom except the ones of the metric. This means the quintessence models or other entities 
of similar origin are completely out of the consideration presented below. 

Under the conditions formulated above, the quantum contributions can be taken in 
form of a power series in the derivatives of the metric. In the cosmological setting, the 
second order in derivatives term means there are only the 0(H 2 ) - like contributions 2 . 
According to the consideration presented in [39], this means that the quantum corrected 
vacuum energy density p\ = A/(8irG) and the Newton constant satisfy the equations 

PA = Co + 7— Mp p 2 , v 

( P + PA )dG + Gdp A = 0, 
3H 2 

P + PA= 8^G> 

2 Let us note that linear in H noncovariant terms are not favored also by the analysis of observational 
data [49]. 



a M 



12 7T M 2 



(2.1) 
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where p is the energy density of matter, fi = H, while v and M are some undefined 
parameters dependent on the unknown quantum corretions. 

The solution for G = G(H; v) can be easily found to be (see also discussion in [50]) 

where G(Hq) = Gq = 1/Mp is the initial value of G. It is interesting to note that the 
relation (2.2) implies that the Appelquist - Carazzone decoupling does not apply to the 
Newton constant. The possible reasons for that has been explained in [39] and we will not 
repeat them here. 

Let us look for a natural extension of the relation (2.2). The presence of an arbitrary 
dimensional parameter \x in Eq. (2.1) makes the formula (2.2) much more general than the 
purely cosmological setting [39] . Indeed, both sorts of hypothetic quantum contributions 
come from the Feynman diagrams which have some number of external legs of a back- 
ground metric. The internal lines of these diagrams may be of the matter fields in case 
of semiclassical approach or also include quantum gravity propagators and vertices in case 
we deal with some model of quantum gravity. 

In all cases the external legs of metric are characterized by some dimensional parameter, 
say the typical energy of the gravitational field in a given physical situation. This energy 
should be associated with the dimensional parameter [i which we have introduced before. 
In the cosmological setting, as far as the process of our interest is the expansion of the 
universe, the typical energy of the gravitational field can be associated with the Hubble 
parameter, so it is natural to set /i = H [46, 39]. However, in another physical situation the 
typical energy parameter \i may have another natural identification. From this perspective 
the relation (2.2) means we can expect similar relation in the general case, where the Hubble 
parameter will be replaced by some, yet unknown, combination of the metric components 
and their derivatives. In the present work we are interested in the application of the general 
relation 

G ^ = n w 2/ 2^ > ( 2 - 3 ) 

l + i/ In (fJ, 2 /^) 

to the rotation curves, and therefore are interested to use the phenomenologically best 
and also natural choice for \i in the corresponding setting. In the previous papers on the 
subject [41, 53, 39, 51] the choice of the scale was done according to [i ~ 1/r, where r is 
the radius of rotation for a given star. It is obvious that this choice of scale is not really 
natural from the quantum field theory viewpoint, because the inverse distance can not be 
seen as a parameter which characterizes the energy of the gravitational field quanta in the 
external leg of the corresponding Feynman diagrams. The inverse distance is a parameter 
which does not depend directly on the intensity of the gravitational field. 

From the observational viewpoint, the above choice of scale leads to an additive con- 
stant contribution to the rotation velocity of galaxies [41, 53, 39, 51]. It is dubious that 
an additive constant to the circular velocity, which moreover is the same for every galaxy, 
can remove or alleviate the need of dark matter on the modeling of galaxy rotation curves 
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(see our forthcoming comments, and [69] for a criticism on this approach regarding another 
model with similar phenomenological consequences). On the other hand, it is remarkable 
that quantum corrections may lead to a tiny variation of G which, in turn, may contribute 
to galaxy rotation velocity in an observable way. 

From the consideration presented above it is clear that a natural parameter to be 
associated with \x should be something which is directly related to the intensity of the 
gravitational field. Indeed, we can not know this intensity beforehand, since it depends on 
the Newton "constant" , which is supposed to vary. In the present paper, in the description 
of the gravitational field at the galaxy scale, we shall choose the Newton gravitational 
potential as the relevant parameter and set 



where &o and a are in this paper presented as phenomenological parameters and <3? Nc „t is 
the Newtonian potential computed with the boundary condition of it being zero at infinity. 
Due to the logarithm in the G(fi) expression, the physical consequences of the above setting 
are essentially insensitive to the value of $o> as long as its value is not several orders of 
magnitude different from a typical value of $ Ncwt in the system under consideration, which 
we will always assume. The dimensionless parameter a has an important role, it must 
be correlated to the mass of the system under consideration. From the rotation velocity 
formula for this model, to be seen in eq. (2.18), minimum consistence requirements show 
that a must go to zero when the mass of the system under consideration goes to zero. 
Considering the application to galaxy rotation curves, the growth of a with the mass is 
also expected. The ratio V^ ewt /$ Newt is independent of the total (baryonic) mass of the 
system, while the (baryonic) Tully-Fisher law [52, 2] states that the observed rotation 
velocity at the outer parts of disk galaxies (i.e., their non-Newtonian part) increases with 
the galactic mass. From eq. (2.18), the latter statements imply that a must increase with 
the galactic mass. We hope that future developments of the theory here presented will 
unveil the structure of a, or justify it as a reasonable effective parameter. 

Contrary to Newtonian gravity, the value of the Newtonian potential at a given point 
does play a significant role in this approach. This sounds odd from the perspective of 
usual Newtonian gravity, but this is not so from the General Relativity viewpoint, since 
the latter has no free zero point of energy. In particular, the Schwarzschild solution is not 
invariant under a constant shift of the potential. The last argument shows that in making 
the choice (2.4) we are implicitly using relativistic arguments. In what follows the reader 
will see that the relativistic effects really play an essential role in our model. 

Solar system and laboratory experiments on gravitation have not yet reported any de- 
viations from (classical) General Relativity. The above setting is also consistent with these 
observations at smaller scales if a decreases sufficiently from its value for the Galaxy to its 
corresponding value for the Solar system 3 . This argument only shows the existence of a way 

3 Since the baryonic mass in our Galaxy is about 10 orders of magnitude greater than the Solar system 
mass, it is reasonable to assume that for the aforementioned experiments one would find /i//io ~ 1, and 
hence G(fi) ^ Go- 




(2.4) 



G 



for finding compatibility with the known tests of General Relativity, but the final answer 
on the form of quantum correction can be achieved only through the direct calculation of 
these corrections from one side and the detailed investigation of their phenomenological 
consequences on another one. 

2.2 The dynamics for a slowly varying gravitational coupling parameter 

Let us consider the dynamics of particles on a theory with weak variation of the gravi- 
tational coupling parameter G. We will keep in mind the relation (2.3) and pay special 
attention to the scale identification (2.4), but until some point will write formulas in a most 
general form, for G close to Go and weakly dependent on the space coordinates x l , namely 
GdGjdx 1 ~ GqOG / dx l . We call it a slowly varying gravitational coupling parameter 4 . 
Our starting point will be the Einstein-Hilbert action 5 

s = c 7 d4 ^T^r < 2 ' 5 > 

where g is the determinant of the metric, R is the Ricci scalar, and G is the variable 
gravitational coupling. Since the cosmological constant A does not play an essential role 
in the rotation curve analysis, we simply set it to zero. 

In what follows we will need a nonrelativistic approximation to the Einstein equations 
and, mainly, to the equation of motion for the test particle in a theory with variable G. 
First of all, let us restrict the consideration by a weakly varying G, such that 

G = G + 5G = G {1 + k), |k|«1. (2.6) 

In the last equation Go is the nonperturbed constant value of the Newton gravitational 
coupling. In particular, we will see in the next sections that the appropriate value of the 
parameter v in Eq. (2.3) is about 10 -7 . This is the value which provides a good fit of the 
rotation curves without the addition of dark matter and, definitely, it is very small. As a 
consequence we can always keep only the first order in the expansion into series in v. The 
same will be done with the general ratio k = 5G/Gq. 

In order to obtain the expression for the corrected gravitational potential in the non- 
relativistic approximation, we can do a straightforward variation of the action in respect to 
the metric, or simply follow the method of [39]. In order to link the metric in the variable 
G case with the standard one, we perform a conformal transformation in the action (2.5), 
according to 

G 

9iiu = -q 9iw = (1 - k)9vu, ■ (2.7) 

4 From the G(/x) expression, one can check that it is indeed a slowly varying function in the mathematical 
sense, but the main reason for the use of this name comes from the approximations G(n) is expected to 
satisfy. 

5 We use ( — H — h) as the space-time signature. For clarity, considering the application to rotation curves, 
the velocity of light c will always be explicit. 
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It is easy to see that, up to the higher orders in k, the metric g^ v satisfies the usual Einstein 
equations with constant Gq. The nonrelativistic limit of the two metrics g^ u and g^ is, 
therefore, 

.2* 2$ Nowt 

goo = -1 2* and 9oo = -1 9— , ( 2 -8) 

c c 

where $ Ncwt is the usual Newton potential and $ is yet unknown potential corresponding 
to the solution of the modifies gravitational theory (2.5). It is easy to see that, in view of 
(2.6), we have 

300 = -1 - ~g = (1 + «)c/00 = (1 + ^—) ~ -1 ^ K ( 2 - 9 ) 



and, hence, 



c 2 c 2 5G 

$ = $Ncwt + 7T « = $Ncwt + -iTPT- • (2.10) 
2 z Ctq 



For the nonrelativistic limit of the modified gravitational force we obtain, therefore, 

c 2 G' j 

= -< wt - — , (2.11) 
2 Go 

where we used the relation G' 1 = (5G)' 1 . In the last formulas and in what follows, all indices 
after a comma denote partial derivatives, while semicolons denote covariant derivatives. 

Another way to arrive at the same formulas is as follows. The affine connections which 
correspond to the two metrics are related as 

f£ A = Kx-^~ o {%G,x + KG,, - g vX G>») + 0(k 2 ) , (2.12) 



where 

T^a = 2^ ^9ua,\ + 9\a,u ~ 9uX,a) and = -g^ (g ua , x + g\ a , u ~ 9v\, a ) 

Since a particle of mass m couples to gravity through the action 



mc / ^-g^dx^dx", (2.13) 



the gravitational acceleration felt by a nonrelativistic test particle in the static weak field 
limit is 

£ = - r -^f w - r » = - f » - k (-« g; ') • <2 - 14) 

where r is the proper time, dr 2 = —g^dx^dx" , with x° = ct. It is easy to see that the last 
equation is equivalent to (2.10). 

Finally, replacing dr 2 = c 2 dt 2 in (2.14), we arrive at the nonrelativistic equation for 
the gravitational acceleration of a test particle in the form equivalent to (2.11), 
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Let us note that the Newtonian potential in the equation above is not exactly the 
Newtonian potential found from Too of the original theory. However, it is easy to see that 
the difference is very small for a very small value of k. The relevant energy-momentum 
tensor is that of a fluid of negligible pressure, namely 

ydust dust t r j j 

1 flU ~ P U fl U U > 

where = dx^/dr is the four- velocity with dr 2 = —g liv dx il dx v '. Furthermore we know 
that p dnst scales with the conformal transformation as oc (Gq/G) 3 ^ 2 , at least up to the 
0(ft)-terms. Hence, 

G_ 



Q~-/^ Tidust O^/^i / rpdust Q^S~1 rridust 

^0 



This approximation is a very good one, since the observational uncertainties on p dust are 
much greater than (k/2) p dust . In conclusion, within a very good approximation, the New- 
tonian potential that appears in eq.(2.15) is the one induced by the energy density of the 
original theory. 

The immediate consequence of the acceleration formula (2.15) for an axisymmetric 
system with radial coordinate R is that the rotation velocity V is given by 

V 2 *V 2 cvt + R C ^. (2.16) 
I Go 

The dimensionless ratio ^ gives the fraction of c 2 which should be added to the square 
of the Newtonian circular velocity. Hence, even the tiny variations of the gravitational 
coupling G across a galaxy are sufficient to generate significative changes in the observed 
rotation curve. 

Fixing the variation of G as in (2.3), the relation (2.16) becomes 

V 2 = V 2 e ^-vc 2 R^. (2.17) 
A* 

Before starting the analysis of the rotation curves on the basis of eq. (2.17), let us 
make one more observation concerning the previous attempts on the fi scale setting. From 
eq.(2.17), the effect of setting \i oc R~ a , as suggested in [39, 51, 53] (in these references 
with a = 1), is to add a constant contribution proportional to a to the circular velocity. 
Considering that this new effect has a significant contribution to the rotation curves of 
galaxies at the farthest observable regions from the galaxy center (large R), the least one 
should expect from this constant is that it should vary from galaxy to galaxy (in order to 
replace dark matter). Otherwise, there is no hope to explain the diversity of rotation curve 
velocities at large R (e.g., see the sample of galaxies studied in this paper). In conclusion, 
following this approach, a must vary from galaxy to galaxy. The problem with the setting 
/i oc R~ a , when applied to real galaxies, appears in the rotation curve behavior close to 
the galactic center. At this region, there are many galaxies that display a rising curve 
that starts at few km/s. Hence, if a constant velocity is added to help the large R part of 
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the rotation curve, the same constant will be troublesome to the central part of the same 
galaxy. 

Another approach to the setting fi oc Rr a comes from a direct integration of the 
effective potential of a point particle, instead of using eq. (2.17). First of all, this procedure 
of summing the potentials constitutes an extra assumption. Second, the integrand has a 
term that, at most, only depends on the mass trough a, thus the meaning of this integration 
should be analyzed with care. If a has the same value for all the points that constitute the 
galaxy, then the circular velocity will depend on a new parameter, the model cutoff for very 
large R, which we label Rf. In this case, the resulting rotational velocity is essentially the 
Newtonian one added by a term proportional to R 2 /Rf- Considering the innermost region 
of a galaxy, this R 2 behavior is much better than the previous one in which a constant 
is added to the Newtonian rotation curve. On the other hand, it is hard to match this 
quickly increasing behavior to the outer parts of many galaxies. At last, if a depends on the 
local matter density, its mass dependence for each point should be close to a square root to 
account for the Tully-Fischer law when galaxies are seen far away (similarly to the approach 
MOND uses), but this non-linear dependence on the mass is not in conformity with the 
summation of potentials; therefore this approach also fails. In conclusion, // oc R~ a is 
not a satisfactory setting for real galaxies, both from a theoretical perspective and from a 
phenomenological one. 

Finally, we can conclude that eq. (2.4) represents a more natural dependence between 
the gravitational coupling and the gravitational potential in the weak field limit that can 
account for the missing matter problem in galaxies. 

Using the setting of fi as given by eq. (2.4) and = vac 2 , the expression for the 
circular velocity becomes simply 



where the relation V 2 evrt = R($ Ns wt),R was used. 

For a point particle, V^ ewt /$ Newt = —1, hence the above velocity profile describes a 
rotation curve around a massive point particle that is Newtonian in its innermost region 
and becomes flat with velocity Vqo in its outermost region. If matter is the source of gravi- 
tational force, than Voo must have a (direct or effective) mass dependence, and in particular 
it must go to zero when the mass goes to zero. Such dependence is also expected from 
the Tully-Fisher law, since it roughly states that the square of the galaxy rotation curve 
velocity scales with the square root of the galaxy luminosity (which in turn is proportional 
to the galaxy stellar mass). 

One may use the Tully-Fisher law, or variations (e.g. [2]), to further constrain the 
model, possibly unveiling a suitable dependence of on other galaxy parameters (e.g., 
mass, morphology, age...). However, contrary to the original formulation of MOND, we 
will let V^o to be a free parameter which may in the future be fixed. 




(2.18) 
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3. Galaxy rotation curves 



3.1 General properties 

The stellar component of many disk galaxies can be fairly described by single thin disk 
whose matter density decays exponentially with the radius [54], 

where M* stands for the total stellar mass and Rd is the disk scale length. Usually, a 
significative improvement can be found by using a thick disk with matter density given by 

p ^R }Z ) = a ,(R) Sech 2 ( - ) , (3.2) 

with zq = Rd/5 [55, 11]. In the absence of a relevant gas contribution, from the above 
matter distributions one infers the corresponding Newtonian gravitational potential and 
hence the velocity profile of the model (2.18). 

In this subsection the non-Newtonian contribution to the square of the Newtonian 
velocity, as given by eq. (2.18), 

Vl R (ln|$ Ncwt |) ii? , (3.3) 

will be compared to the Isothermal profile contribution. The latter is one of the phenomeno- 
logical dark matter velocity profiles that better describes the observed rotation curves (e.g. 
[10, 11, 19]), it was studied in many works, and it also has a non-null asymptotic velocity 
for large galactic radius. The latter property, which is not satisfied by all dark matter 
profiles in agreement with observations (e.g., [16]), makes the Isothermal profile a natural 
one to compare with the model here proposed. 

The velocity rotation curve of a galaxy modeled with the aid of the Isothermal profile 
is given by V 2 = 1/ N 2 cwt + V 2 a , where 



Vl(R) 
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(3.4) 



and V 2 o00 = 4 ir Go po R 2 C . This profile is commonly parametrized by R c and po; but, in 
this section, since the model previously introduced also has a non-null asymptotic velocity 
V^oo , it is convenient to parametrize it by R c and V laoOQ . The Isothermal profile has no direct 
dependence on the baryonic matter distribution, on the other hand a number of scaling 
features that depend on the baryonic matter properties appears when this profile is applied 
to real galaxies (e.g. [1, 18, 17, 13]). 

In Fig.(l) we plot the dimensionless ratio of the additional contribution (3.3) by V^, 
where the Newtonian potential is computed from both a thin and a thick stellar exponential 
disks. In the same figure we compare the results with the Isothermal profile contribution, 
also divided by V^, for some values of its parameters. Considering a galaxy with negligible 
gas mass and whose stellar component can be approximated by an exponential disk, Fig.(l) 
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reveals some properties of our model, including: i) the presence of an effective core of 
radius 6 R ~ Rr>', ii) the existence of a rotation curve maximum at R ~ 5Rd, leading to the 
an approach of the asymptotic velocity from above (contrary to the Isothermal profile); Hi) 
the Isothermal profile can be set to be practically identical to the slowly varying G model 
up to about 5Rd] iv) the thickness of the stellar disk slightly modifies the stellar rotation 
curve, which in turn also slightly modifies the non-Newtonian contribution, generating 
a bigger effective core, and displacing the maximum of the non-Newtonian contribution 

Regarding the item i, the correlation be- 
tween the core radius and the disk scale length 
is one of the already known correlations be- 
tween the cored dark matter profiles and the 
baryonic matter [18]. From the item ii above, 
a galaxy with negligible gas mass and expo- 
nential stellar density profile, whose rotation 
curve clearly increases for R > 6Rd, cannot 
be satisfactorily described by the the variable 
gravitational coupling model. To our knowl- 
edge, no regular galaxy with these proper- 
ties have ever been found. Among our sam- 
ple of galaxies, the galaxy NGC 2841 which 
can be seen in Fig. (5) is the one with the 
poorest gas content. It does has a relevant 
bulge component, but its effect is small for 
large radii. Its rotation curve was measured 
by more than 10 disk scale lengths, and it 
clearly displays a decay with the radius in 
agreement with the above properties. 

Regarding the effects of adding a gas component, Fig. (2) presents an example of 
the deformation caused by the gas on the non- Newtonian contribution (3.3). Details of 
the general structure of the gas rotation curve are still under debate (see [71] for some 
possibilities), but we can fairly illustrate such effect by using a re-escaled version of a 
"typical" gas rotation curve, say the one from NGC 3198 [11]. To this end it is convenient 
to work with dimensionless quantities, thus we use the following dimensionless Newtonian 
potential 

MR/Rp) f ^{R/Rp) 

M*G J Af gas G ' 
where $ gas and M gas are deduced from the NGC 3198 data, while / weights the gas influence 
to the potential. In particular, / = 1 represents a galaxy whose stellar mass is equal to its 
gaseous mass. Both theoretically and from the numerical result of Fig. (2), one sees that 
the greater is the gas fraction /, the smaller is the non-Newtonian contribution to the inner 
6 The latter is determined by finding the Isothermal core radius R c value such that the V^ inflection 
point happens at the same radius of our model. 
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Figure 1: A comparison between the non- 
Newtonian contribution induced by our model 
to the dark matter rotation curve generated 
by the Isothermal profile. Black lines: the 
non-Newtonian square velocity contribution di- 
vided by for thin (solid line) and thick (dot- 
ted line) exponential stellar disks. Gray lines: 
VfjVl for V£ ooc = 1.6 Vl with R c = 0.95 R D 
(solid line) and R c = 2 Rd (dotted line). The 
internal plot is a zoom in the region R/Rd < 
1.2. 



- 12 - 



parts of the galaxy, and the higher is its contribution to the outer parts of the galaxy. This 
feature is in agreement with observations, since the gas-rich galaxies are the Low Surface 
Brightness (LSB) galaxies, which typically display a slowly increasing rotation curve up to 
the farthest observed point. 

Figure (2) shows that rotation curves modeled with our model are more susceptible to 
the gas irregularities than the rotation curves modeled with dark matter (with the Isother- 
mal, Burkert or NFW profiles, for instance). If the gas rotation curve irregularities are too 
much amplified, this can pose a serious problem to the model, since strong irregularities 
may be inserted in the resulting rotation curve (see, for example, the "HI scaling model" 
[68, 10] case). Currently it is hard to say if our model amplifies such irregularities either 
more or less than MOND, since this behavior varies from galaxy to galaxy, depending on a 
number of factors. The results of our model, to be present in the forthcoming subsection, 
show that such reflection of the irregularities does not prevent our model to accommodate 
considerably well the observational rotation curves. It would be interesting to do a tougher 
investigation on this issue by using a larger sample of galaxies and looking for significative 
systematic effects. 

Figure (2) unveils another property of 
the model. Independently on the value of 
/, all the curves in Fig. (2) pass by a single 
point. From (3.3), one realizes that the non- 
Newtonian contribution at a distance Rq, that 
satisfies (3>*/3> gas ) R (Ro) = 0, is not affected 
by the value of /. The existence and posi- 
tion of such point cannot be proved in gen- 
eral due to the diversity of the gas rotation 
curves. However, in the innermost regions of 
a galaxy, the ratio <E>* / $ gas is expected to de- 
crease considerably with the galaxy radius; 
while, in its outer regions, the gas contribu- 
tion typically starts to dominate. The latter 
leads to an increase in the ratio / $ gas and 
thus, under natural assumptions, it implies 
the existence of the point Rq at the outer 
parts of many galaxies. 

3.2 Mass decomposition: our results on 
real galaxies in comparison with other models 

To test the model here presented, which hereafter we label RGGR, we use recent high qual- 
ity data from nine regular spiral galaxies obtained from two different sources [11, 10] and 
compare, for each galaxy, our results with those of three other models, namely: the Isother- 
mal profile, the Modified Newtonian Dynamics (MOND) [20] and the Scalar- Tensor- Vector 
Gravity (STVG) [56]. The first being a well known and phenomenologically successful 
profile for dark matter; the second is the most well known, historically important and con- 




R/R D 



Figure 2: An illustration of the deforma- 
tion on the non-Newtonian contribution in- 
duced by the gas presence. Computed us- 
ing a re-cscaled version of the NGC 3198 gas 
data from THINGS. The solid black line is 
the squared non-Newtonian velocity contribu- 
tion divided by , for a galaxy with thin stel- 
lar disk and no gas (/ = 0). The gray lines are 
the same as before but with the following values 
of /: 0.2, 0.7, 1.2, 1.7, 9.7. 
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siderably successful modified Newtonian gravity model; and the third is a recent General 
Relativity modification that in its weak field limit has similarities with MOND, while it 
is also capable of providing good results on lenses and other General Relativity related 
phenomena (see [57, 58] for recent reviews). 
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The models 

The free parameters of each model According to our model, galaxy rotation 

curves can be explained without the need of 
dark matter, once renormalization group correc- 
tions to General Relativity are taken into con- 
sideration. Besides dark matter profiles, a natu- 
ral model to compare our results with is MOND, 
since it is currently the most influent and per- 
haps the most successful approach to explain 
the missing matter problem in galaxies without 
dark matter. MOND will be here tested only 
in its original formulation [20] , and possible im- 
provements on its results due to small changes 
in either galaxy inclination or the distance will 
not be considered. The details of such changes 
to MOND can be found in [59, 28]. To allow for 
a comparison of model results using the same 
data and procedures, MOND will be here fitted 
considering the data from Sample A [11] (this 
is the first time that MOND is fitted with this data). MONDian fits using the data of 
Sample B [10] were already done, see [10, 28] for instance. Nevertheless, we will do all 
the procedures of the fits for MOND to this Sample too, in order to have uniformity on 
the results and to allow a fair comparison of the results. We remark, in particular, that 
in the original work on the subject [10] the minimization procedure is different to the one 
employed here (we use the most usual convention for x 2 ; soon to be specified). 

The MONDian acceleration, already written in terms of its circular velocity, is given 

by 



Table 1: (1): "RGGR" labels the model we here 
propose based on Renormalization Group corrections 
to General Relativity, "Iso" labels the dark mat- 
ter Isothermal profile, "MOND" the Modified New- 
tonian Dynamics and "STVG" the Scalar- Tensor- 
Vector Gravity model. (2): Lists the model's pa- 
rameters that are free to vary from galaxy to galaxy. 
(3): Lists all the parameters whose values are de- 
duced from the fitting procedure. Mp and M B label 
the disk and bulge stellar masses respectively. For 
galaxies with no bulge, M B is not a free parameter. 
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R clq I R 

" 8 cm/s 2 [60]. Consequently, given the function 



with n(x) = x/Vl + x 2 and ao = 1.2 x 10 
V Newt (R), V MOm (R) is at once found. 

Likewise the model we present in this paper (RGGR), the STVG model also uses the 
running of the gravitational coupling motivated by a renormalization group approach to 
gravity. Its consequences were studied for a large group of galaxies in 7 [25]. The STVG 



7 Reference [25] centers its analysis on the Metric-Skew- Tensor-Gravity (MSTG) model, which neverthe- 
less induces the same acceleration law for weak gravitational fields (see also [56]). MSTG also uses the 
running of G to explain galaxy dynamics. Since here we only analyze galaxy rotation curves, we could label 
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model is a natural model for us to compare with. To this end, this model will be fitted 
with the same data and procedures we use to fit our proposed model; in particular we will 
use the data from Samples A [11] and B [10] (this is the first time that this model is fitted 
with this data). 

The circular velocity induced by the STVG model is [25] 



where M Total is the total "baryonic" mass of the galaxy given by M* + M gas , M$ and ro 
are related by Mq = Tq aQ TVG /Go, with ao TVG = 6.90 x 10 -8 cm/s 2 . For galaxies whose 
outermost measured velocity lies at a radius greater than 12 kpc, ro is set to be 13.92 kpc; 
otherwise, ro is 6.96 kpc [25]. From the theory, ro is supposed to vary from galaxy to 
galaxy following the renormalization group flow, but these two values are sufficient for a 
good agreement [25]. 

For the Isothermal profile, we use the velocity given in eq. (3.4) parametrized by R c 
and po, as usual. And, for the RGGR model, we use the velocity given by eq. (2.18), but 
parametrized by a, where Voo = a c 2 10~ 7 ; that is, for convenience, v is fixed to be 10~ 7 . 
In the end, the parameter that matters for galaxy rotation curves is the product au. 

Table 1 lists the models and its corresponding free parameters. Besides the parame- 
ters that are intrinsic to each model, the disk and bulge stellar masses also appear as free 
parameters. This is clarified in the next subsection. 

Stellar masses and the fitting procedure 

To derive the theoretical rotation curve for each galaxy and each model, there is a 
number of relevant parameters whose values are not known beforehand. For instance, 
regarding the Isothermal profile, each galaxy can have its own value of core radius R c 
and central density po, while for the RGGR model the parameter a is free. These are 
parameters that are model specific. Moreover, although the stellar mass distribution of 
regular disk galaxies is expected to follow approximately an exponential profile [54], the 
total stellar mass is subject to a considerable uncertainty. The main trend is an increase 
of the stellar mass with the increase of the absolute luminosity of the galaxy, hence we use 
the stellar mass-to- light ratio in place of the stellar mass as a galaxy parameter. One can 
always convert the latter into stellar mass by using the luminosity values that are provided 
in Table 2. 

We use the standard procedure for handling the stellar mass-to-light ratio. In a first 
step, it is considered to be a free parameter whose value can be found in conjunction with 
the selected model by a least-square fit, in which the observational velocity and its error 

this model to be tested as either STVG or MSTG. — See however "Note added". 




(3.6) 
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bars are considered. Namely, a numerical search for the global minimum of 



All observed points . 0ba , . s 2 

e ( ' " E l:r ) (a7 » 

is done for each galaxy and each model. In the above, i runs all the observed points of a 
given galaxy, I^° bs and Error, are respectively the i-th observed rotational velocity at the 
radius Ri and its error. The quality of the rotation curve fit (considering its shape alone) is 
characterized by y^ ed = Xmih/ (Number of observed points — Number of parameters being fitted), 
where xLin * s the minimum value of x 2 found in the previous procedure. This procedure 
yields a value to the stellar mass (and hence to the stellar mass-to-light ratio) and to 
additional model parameters if present. The associated 1 a errors to the fitted parameters 
are deduced by searching for both the maximum and minimum of the same parameter 
inside the region x 2 < xtiin + 1> where all other parameters are free. 

The galactic stellar mass is not a true free parameter, since estimates for its value 
can be derived from stellar population synthesis models, e.g. [61, 62]. After the fitting 
procedure, it is analyzed if the resulting mass-todight ratio is a satisfactory one or not. 
The expected values are listed in Table 2. 



About the selected samples 

In order to analyze the RGGR model, and to compare it to other proposals, we selected 
two recent high quality samples of spiral galaxy data. The data used can be found in[ll] 
(Sample A) and [10] (Sample B). Not all galaxies of these references are here fitted. From 
the sample of galaxies fitted in [11], we select five whose observed rotation curve are most 
regular and that the stellar component can be very well approximated by one or two thick 
exponential disks, each disk with density given by eq. (3.2). Following the conventions of 
[11], if two exponential disks are used, the disk that dominates the inner part of the galaxy 
it is called "bulge". The standard modeling of the bulge with the i? 1 / 4 profile seems less 
suitable for this sample of galaxies [11]. 

The regularity of the observed rotation curves is an important issue. Upon avoiding 
galaxies with either ill understood or just exotic features, emphasis is given to those that 
the presented models (built under assumptions of regularity) have a chance of handling 
and providing meaningful answers. 

The same procedures used in [11] for the fitting we also employ here, with the single 
exception that we use thick exponential disks as approximations for the complete stellar 
density (instead of only large radii continuations from the photometric data). The expo- 
nential disk approximation is a robust one [54]. The details on the stellar rotation curve 
can either improve the concordance with observations (e.g., NGC 2403) or worsen it (e.g., 
NGC 3198). 

For the Sample B we use all the conventions used in [10], apart from their x 2 min- 
imization. From the five galaxies presented in [10], only one (ESO 79-G14) is not fitted 
here since it does not have an exponential stellar profile. 
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Table 2 displays general information on the galaxies that are part either of Sample A 
or B. The corresponding values of Rp and Rb explicitly cited in [11] are not corrected for 
inclination, and, when cited, refer to the exponential continuation only (except for DDO 
154). The values displayed in this table were found from fitting each stellar component 
rotation curve with the rotation curve of a thick exponential disk, in which zq is set to be 

1 /5 of the corresponding exponential disk scale length. The quoted values of the mass-to- 
stellar ratio are only used to compare with the values found from the fitting procedure. 

The Newtonian rotation curves 

The gas contribution to the Newtonian rotation curve can be evaluated from the New- 
tonian potential of the gas, which can be straightforwardly deduced from the gas surface 
density. (There are different ways to compute the Newtonian rotation curve, but since the 
RGGR model depends explicitly on the Newtonian potential, this seems to be the most 
suitable one.) The surface density of the neutral atomic hydrogen (HI) is inferred from 
the 21 cm radiation of each galaxy (see [67] for Sample A and [10] for Sample B), and it 
differs from the gas surface density by a multiplicative constant (below specified). In this 
work, our starting point for the gas rotation curve is the neutral atomic hydrogen surface 
density, as provided by the respective collaborations. 

The adopted procedure to find the Newtonian rotation curves goes as follows for Sample 
A: from the surface density of the neutral hydrogen multiplied by 1.4, to account for the 
presence of helium and metals, we deduce its Newtonian potential with the boundary 
condition of it being null at infinity. The latter is then numerically derived to find the gas 
rotation curve. For each stellar disk, we use the appropriate disk scale length from Table 

2 to find the Newtonian potential of a thick exponential disk, from eq. (3.2). Except for 
an overall multiplicative factor for each disk, which is the total stellar mass, the rotation 
curve is obtained by numerically deriving the potential. 

For Sample B, following the conventions used by the original reference, we use the same 
procedure as above, except that the surface density of the neutral hydrogen is multiplied 
by 1.33, to account for the presence of helium, and that the stellar disks are thin; hence 
the Newtonian potential is found from eq. (3.1). 

Some galaxies, close to their center, have a gaseous rotation curve component with 
a negative contribution to the total circular velocity. Effectively, this leads to V? m < 
for a certain radius interval, which can describe a physical rotation curve at that radius 
if its sum to the other contributions results in a positive total squared circular velocity. 
In the plots, it is customary to symbolically express this property in the following way: 
V ew = if V g l > 0, and V gaa = -ypi%. otherwise. 

The results 

The results of the fitting procedure can be found in Tables 3 and 4, and in the Figs. 
(3, 4, 5, 6, 7, 8, 9, 10, 11, 12). 

Comparing the values of a from tables 3 and 4 to the galaxy luminosities in table 2, 
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one can clearly note a correlation of a with the galaxy luminosity, as already expected from 
Tully-Fisher arguments. Since the galaxy luminosity correlates with the galactic mass, we 
see that a increases with the galactic mass, as foretold. Future analysis of this model with a 
larger sample of galaxies should constrain the dependence of a on the galactic parameters. 

Firstly we warn that the absolute values of the xLi are n °t of considerable physical 
significance, since the most relevant information the error bars carry is their relative sizes. 
There are differences on conventions regarding the error bars values between the Samples 
A and B [11, 10]. The relevance of Xied nes on the comparison of its values between 
different models for the same galaxy. The RGGR model typically scores almost so well 
as the Isothermal profile and considerably better than both MOND and STVG. There are 
however two exceptions, the galaxies NGC 3621 and ESO 116-G12. The RGGR Xr Cd value 
is 3 and 4 times greater than the Isothermal ones for the respective galaxies. For NGC 
3621, RGGR has a xl d value as good as STVG, but worse than MOND; while for ESO 
116-G12, the RGGR model scores worse than STVG and better than MOND. 

The RGGR problem with NGC 3621 lies in the middle radius range, while it can fit 
well both the rising curve and the irregular shaped final part, reflecting the gas behavior. It 
is very far from a disastrous result, and, moreover, the result this fit yielded for the stellar 
mass-to-light ratio is the one closer to the expected Kroupa value. For ESO 116-G12 the 
problem is not localized, but the final gas wiggle is, for this galaxy, a difficulty. 

For the stellar mass-to-light ratios, the model that systematically provides the best 
values, according to the expected values quoted in Table 2, is the RGGR model. 

For the Sample A, the RGGR model yielded twice a mass-to-light ratio that is outside 
the 50% range of the expected values. One for DDO 154, a galaxy that no model can yield 
a good mass-to-light ratio, but that the RGGR model yields the closest one. And another 
for the bulge of NGC 2841, which is slightly above the allowed range. Probably a very 
small decrease in the Rb value could solve the previous issue and also improve the rotation 
curve concordance for small radii. 

On comparing our Isothermal results with those of [11] care should be taken because 
we use exponential approximations for the complete stellar components, an approximation 
that [11] only uses for DDO 154. However, our results are considerably alike. The small 
structures of the stellar rotation curve can either improve the fit (e.g., NGC 2403), or 
worsen it (e.g., NGC 3198). 

In the Sample B only MOND and STVG found mass-to-light ratios outside the ex- 
pected range, while the Isothermal was three times on the verge of the expected values. 
The shapes of the curves values found for the Isothermal and MOND models are essen- 
tially the same of those found in [10], but there are systematic differences on the inferred 
values of the fitted parameters, which are expected since contrary to [10] here we used the 
standard (or the most common) x 2 definition for the minimization procedure. 

The Isothermal profile has parameters strong correlated to be meaningfully fitted in 
some cases. Two considerably different set of parameters can yield the same x 2 B d within 
two precision digits. As an extremum example, we cite an issue we found on NGC 2841. 
The result that we presented in Table 3 has x 2 = 17.9, which within three precision digits 
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is not truly the global minima of x 2 ■, since a completely different set of parameters resulted 
into x 2 = 17.7. The parameters of the latter fit are Y*d = 1.2 M /L , Y*b = 1-1 M & /L Q , 
R c = 7.3 kpc and po = 25 x 10~ 3 Mq/pc 3 . Note that both stellar mass-to-light ratios have 
changed from good to far from the expected values. The result presented in Table 3 on 
NGC 2841 was the only one we selected considering both the x 2 minimization and the 
likeness of the mass-to-light ratio. 

We did some experiences by letting ao in MOND to be a free parameter. The fittings, 
considering the shape of the curve, can be significantly enhanced in this picture, but no 
significant improvement on the mass-to- light ratios was found. For the latter, it is necessary 
to change the shape of fJ,(x), likewise done in [28]. 

4. Conclusions and discussions 

We presented a model, motivated by quantum corrections to the Einstein-Hilbert action, 
that introduces small inhomogeneities in the gravitational coupling across a galaxy (of 
about 1 part in 10 ) and can generate galaxy rotation curves in agreement with the obser- 
vational data, without the introduction of dark matter as a new kind of matter. Considering 
the samples of galaxies evaluated in this paper, the quality of the rotation curve fits and the 
associated mass-to-light ratios from our model is a bit lower, or about the same, than the 
Isothermal profile quality, but with one less free parameter. We also compared the results 
of our model with MOND [20] and STVG[56, 25], and at face value our model yielded 
clearly better results. 

Our results can be seen as a next step compared to the previous models motivated by 
renormalization group effects in gravity, [39, 51]. Their original analyses could only yield 
a rough estimate on the galaxy rotation curves, since they were restricted to modeling a 
galaxy as a single point. Trying to extend this approach to real galaxies, we have shown 
that the phenomenologically optimized proper scale for the renormalization group phe- 
nomenology is not of a geometric type, like the inverse of the distance, but is proportional 
to the Newtonian potential with null boundary condition at infinity. This setting is shown 
to be in agreement with both theoretical expectations and with observations. 

The essential feature for the rotation curve fittings is the formula (2.18), which is 
by itself a remarkably simple formula that provides a very efficient description of galaxy 
rotation curves. We have shown that it can be derived from the assumption of a grav- 
itational coupling parameter with a very small departure from the constant Go, such 
that, in the weak field limit, it depends on the logarithm of the Newtonian potential, 
e.g. G(<I> N ewt) = Go [1 — /3 In ($ Newt /$o)]- I n the last equation, (3 is a positive and very 
small (~ 10~ 7 ) effective parameter that necessarily depends on the distribution of mass of 
the system. The precise relation between f3 (or olv using the theory of quantum corrections) 
with the galaxy parameters is yet to be unveiled. To this end, from a phenomenological 
approach, a larger sample of galaxies would be necessary to find meaningful results. On 
the other hand, using the Tully-Fisher law, it is not hard to guess that /3, for disk galaxies, 
should scale with the mass approximately as \AM, in which M is the baryonic mass of the 
galaxy. The results for a in tables 3 and 4, clearly show that it indeed increases with the 
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galactic mass. In order to avoid premature statements, we think that the nature of f3 (or 
a) should be explored in more detail in the forthcoming papers. 

In the present paper we have addressed the issue of generating galaxies rotation curves 
by developing and extending the proposal of [39]. It turns out that, at least for the 
sample galaxies we dealt with here, one can explain these curves without invoking the dark 
matter concept. Does it mean that we can really be free of dark matter in constructing a 
realistic cosmological scenario? The answer to this question is probably negative. It is well 
known that there many other issues associated with dark matter, or the lack of it, that are 
necessary to take into account. In particular, the results on the density perturbations and 
related issues such as cosmic microwave background radiation (CMBR) and the large scale 
structure (LSS) data, baryon acoustic oscillations, big bang nucleosynthesis, gravitational 
lensing and others (see, e.g. [72] for a recent review). In all these observational and 
experimental issues the standard way of explaining the date is to assume the existence 
of dark matter. At the same time, one can not underestimate the fact that the rotation 
curves may be explained without dark matter effect. One can imagine, for instance, the 
scenario with essentially smaller amount of dark matter which has slightly different set of 
properties compared to the usual CDM model. In this case the rotation curves will be 
explained by summing up the effect of quantum corrections and the one of the dark matter 
content. The preliminary analysis shows this possibility can not be ruled out [73]. Finally, 
there is a chance to address all mentioned issues trading a large amount of cold dark matter 
content by another one, which can have qualitatively other origin, and possibly invoking 
the quantum corrections to the cosmological perturbations spectrum (see, e.g., [64, 65]). 

Note added. After the first version of this work was finished, we became aware, through 
[58], and contrary to the expectations in [56], that there are small and perhaps sensible 
differences between the fits of the rotation curves using either the STVG or the MSTG 
formulations. If precision on the nomenclature and results is at stake, our results on 
the model developed by Moffat and collaborators apply to the older formulation, namely 
MSTG. 
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NGC 3198 
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0.5-1.8 




NGC 7339 
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1.53 




0.5-1.8 





Table 2: (2): Distance of the galaxy (Mpc). (3): Total luminosity (10 10 Lq), the band depends on the sample: 3.6 
fim for Sample A and / band for Sample B. (4): "Bulge" luminosity (lO 10 !/©). (5): Stellar disk scale length (kpc). 
(6): The "bulge" scale length (kpc). (7): Expected stellar mass-to-light ratio for the disk. The values quoted for 
the Sample A arc the Kroupa ones. (8): The same as before, but for the " bulge". We stress here that the expected 
mass-to-light ratios are not used as an input for the fitting procedure, all the fits are done with free stellar masses. 
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Table 3: (1) Name of the galaxy. (2) Model used to do the mass decomposition. (3) Reduced x 2 value. (4) 
Inferred stellar disk mass-to-light ratios in the 3.6fim band (Lq/Mq). (5) The same of the previous item but for the 
"bulge" (Lq/Mq). (6) Value of the dimensionless RGGR constant (V£, = a c 2 10~ 7 ). (7) The Isothermal profile 
core radius (kpc). (8) The Isothermal profile density (10 _3 Mq/pc 3 ). A "t" signs a mass-to-light ratio far from its 
expected value: either 50% greater or smaller than the quoted Kroupa IMF values (considering the fitting errors). 
(*) There is a significant degeneracy on the parameters of this model for this galaxy that is not reflected in the formal 
1 a errors above, see "The results" . 
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Table 4: (1) Name of the galaxy. (2) Model used to do the mass decomposition. (3) Reduced x 2 value. (4) 
Inferred stellar disk mass-to-light ratios in the 3.6/i band for Sample A and in the / band for Sample B (Lq/Mq). 
(5) Value of the dimensionlcss RGGR constant, V£, = a c 2 10~ 7 . (6) The Isothermal profile core radius (kpc). (7) 
The Isothermal profile density (10" 3 M Q /pc 3 ). "f" signs a mass-to-light ratio outside from its expected value (i.e., 
Y* outside the range 0.5 — 1.8 Lq/Mq, considering the 1 a error values). 
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Figure 3: DDO 154 rotation curve fits. The red dots and its error bars are the rotation curve 
observational data, the gray ones close to the abscissa are the residues of the fit. The solid black 
line for each model is its best fit rotation curve, the dashed yellow curves are the stellar rotation 
curves from the bulge (if present) and disk components, the dotted purple curve is the gas rotation 
curve, and the dot-dashed green curve is the resulting Newtonian, with no dark matter, rotation 
curve. 
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Figure 4: NGC 2403 rotation curve fits. See Fig. (3) for details. 
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Figure 5: NGC 2841 rotation curve fits. See Fig. (3) for details. 
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Figure 6: NGC 3198 rotation curve fits. See Fig. (3) for details. 
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Figure 7: NGC 3621 rotation curve fits. See Fig. (3) for details. 
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Figure 8: ESO 116-G12 rotation curve fits. See Fig. (3) for details. 
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Figure 9: ESO 287-G13 rotation curve fits. See Fig. (3) for details. 
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Figure 10: NGC 1090 rotation curve fits. See Fig. (3) for details. 
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Figure 11: NGC 7339 rotation curve fits. See Fig. (3) for details. 
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Figure 12: Stellar disk mass-to-light ratio (Y*d) in the 3.6^m band as a function of the color 
,/ — K, for Sample A. The gray curve depicts the F*£> vs. J — K relation according to the Kroupa 
initial mass function [62, 61] (see also [11]). The J — K color of each galaxy is determined from 
the 2MASS Large Galaxy Atlas [66], except for DDO 154 (see notes in [11]). Each galactic disk 
is represented above by an open circle, with a reference error bar of 50% of the Y*d value (a bit 
bigger than the maximum Y*r> given by the dict-Salpcter initial mass function [61]). The small 
scatter of the circles around the theoretical line is due to small color gradients present in the stellar 
disks. The galactic disks depicted above are from, in order of increasing color: DDO 154, NGC 
2403, NGC 3621, NGC 3198, NGC 2841. The black open squares display the Y* D values and their 
associated errors for each galaxy as inferred from the rotation curve fits for each model. 
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